Quantum computational renormalization in the Haldane phase 
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Single-spin measurements on the ground state of an interacting spin lattice can be used to per- 
form a quantum computation. We show how such measurements can mimic renormalization group 
transformations and remove the short-ranged variations of the state that can reduce the fidelity of 
a computation. This suggests that the quantum computational ability of a spin lattice could be a 
robust property of a quantum phase. We illustrate our idea with the ground state of a rotationally- 
invariant spin-1 chain, which can serve as a quantum computational wire not only at the AfHeck- 
Kennedy-Lieb-Tasaki point, but within the Haldane phase. 



PACS numbers: 03.67.Lx,75.10.Kt,64.60.ae 

Measurement-based quantum computation (MQC) 
proceeds by performing a sequence of single-spin (local) 
measurements on an entangled resource state of a lattice 
or graph. The canonical example of such a resource is the 
cluster state [1] on a 2D square lattice, although recently 
alternatives have been proposed [2-6] . Ideally, such a re- 
source would be natural, appearing as the stable ground 
state of a realistic (experimentally accessible) spin lattice. 
It would also be robust, insensitive to variations in the 
parameters of the Hamiltonian, such that its quantum 
computational ability is attributed to a quantum phase 
in a similar manner to superconductivity and quantum 
magnetism. Evidence of such a quantum computational 
phase has been suggested in a handful of artificial mod- 
els [5-7]. A central problem in this approach, however, 
is that short-ranged variations in a phase, irrelevant to 
the low-energy physics, will in general be extremely dele- 
terious for MQC, where the effect of every single-spin 
measurement is significant for the computation. 

In this Letter, we show how local measurements within 
MQC can transform a ground state in such a way as 
to physically implement a renormalization group (RG) 
transformation, identifying such a quantum computa- 
tional phase and correcting for the short-ranged varia- 
tions. As a specific example, we consider rotationally- 
and translationally-invariant spin-1 chains, which possess 
a Haldane phase containing the Affleck-Kennedy-Lieb- 
Tasaki (AKLT) spin-1 ground state [8]. As shown in [ ], 
measurements on the AKLT state can simulate an arbi- 
trary single-qubit gate sequence in the quantum circuit 
model, i.e. it is a quantum computational wire [4, 10], 
and forms a basic constituent of MQC when such chains 
can be coupled. We show that ground states within 
this phase can also function as quantum computational 
wires by appropriately modifying the AKLT measure- 
ment sequences. This modification can be interpreted 
as quantum computationally simulating a renormaliza- 
tion group transformation, distilling out the long-range 



degrees of freedom which are common to the entire Hal- 
dane phase as shown by various classical algorithmic RG 
methods [11-14]. Although specific to this spin-1 model, 
our result suggests that a similar technique may be ap- 
plicable in any phase for which a known resource state 
is a fixed point of an RG flow. Unlike state filtering 
techniques to distill resource states via local measure- 
ment [3, 6, 10, 15], which are strongly dependent on the 
precise description of the initial state, our method im- 
plements a parameter-independent RG that functions ro- 
bustly within the phase. 

The Haldane phase and logic gates. — The AKLT model 
was originally proposed to analyze the so-called ID Hal- 
dane phase of a spin-1 chain, which displays several char- 
acteristic features (see e.g. [16]) such as a spectral gap 
independent of the system size, a diluted antiferromag- 
netic order often measured by the string order parame- 
ter, and an effective spin-^ degree of freedom (the edge 
state) appearing on the boundary of the chain. A generic 
translationally- and rotationally-invariant Hamiltonian 
with nearest-neighbor two-body interactions on a spin- 
1 chain, which takes the form 

H[fi) = JJ2 ^ ' Sj+i ~ /3(S, • S J+1 ) 2 ] , (1) 
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has a gapped Haldane phase for J > and — 1 < /3 < 1. 
At the AKLT point (3 = — |, each term in Eq. (1) is the 
projector onto the spin-2 subspace of neighboring sites, 
modulo an additive constant and scale factor. Moreover, 
H(— j) is frustration- free, and its ground states have ef- 
ficient matrix product state (MPS) descriptions [17]. 

The ground states of finite chains are nearly fourfold 
degenerate, corresponding to a tensor product of two 
two-dimensional edge states. Each edge state can thus 
be thought of as a qubit. To associate one qubit to a 
chain, say on the right, we may fix the left edge by as- 
suming an additional spin-^ particle which terminates 
the chain, along with an s • S coupling. The encoded 
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Pauli operators then take the form of string operators 
S fe = <J k ® e™ Sk <g> • • • ® e i7rSk , for k e {i, y, z}. This is 
depicted in Fig. 1, which also details the protocol of [ | 
to perform logical operations on the encoded qubit by 
subjecting the ground state to a sequence of single-spin 
measurements. During the MQC, we do not consider dy- 
namics under H(/3); it only specifies a family of resource 
states. However, controlling H(f3) at the boundary re- 
alizes a quantum computational wire without the RG 
methods presented here [18], and may be useful for pro- 
viding protection against errors during computation [9]. 
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FIG. 1: (Color online.) Single qubit operations on a Haldane- 
phase spin chain. A chain of spin-1 particles terminated on 
the left by a spin-| particle encodes one qubit on its right 
edge. Single spin measurements, shown in (a), implement sin- 
gle qubit operations on the AKLT chain, with measurement 
in the basis {\Sg i — 0}}f =1 for some Cartesian basis {ei}^ =1 
leading to n rotation around the outcome axis [9]. Each out- 
come occurs with probability 1 / , a. Fixing a "standard" basis 
{x,y,z} (called {|a;),|y),|z)} in the spin-1 state space), the 
first two outcomes of measurement in a basis {x',y',z} ro- 
tated by 6 around the z axis (spin-1 states {\9) ,\8 + n) , \ z)}, 
for \6) = \ [(l+e" 1 ") |a;) + (l-e~ ie ) \y}]) result in the same 
rotation R z (8) of the qubit, followed by a corresponding 
byproduct n rotation of it around x or y (which can be later 
corrected). The third outcome is just a byproduct z rotation. 
Induced rotations become noisy for j3 7^ — |, but can be im- 
proved by buffering, depicted in (b). Here the left and right 
spins are measured first, and the rotation measurement (mid- 
dle) is only attempted when these are both z. Failing this, 
the middle spin is measured in the standard basis, and the 
attempt is repeated in the next block. Concatenating block-3 
buffering is equivalent to block-9 buffering, as shown in (c). 

As the protocol described in Fig. 1 works for any Carte- 
sian basis, it is reasonable to apply it to the ground state 
of any rotationally-invariant Hamiltonian within the Hal- 
dane phase. Indeed, measuring each spin in the same 
fixed basis implements the resulting 7r rotations, i.e. log- 
ical identity operations modulo Pauli byproducts, with 
unit fidelity at any point in the Haldane phase. How- 
ever, as shown in Fig. 2 (blocksize 1), the gate fidelity 
decays away from the AKLT point when measuring in 
more than one basis. 

Such behavior is not unexpected, owing to the relation 
between rotation fidelity and string operator expectation 
values as detailed in the caption to Fig. 2. R z (8) fidelity 
is related to expectations of S^/ ® cr x and Y> z ® o z (for 
doubly-terminated chains), where x' = xcosO + ysin#, 
and for 9 ^ 2im the former is not guaranteed to take on 



nonzero values throughout the phase. This formulation 
using string operators also connects the perfect fidelity 
of standard basis operations to the infinite localizable 
entanglement length of any ground state in the phase [19] . 
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FIG. 2: (a) Fidelity of 5 rotation, the worst case, and (b) 
buffering success probability versus f3 for blocklengths L — 1 
(no buffering), 3, and 9 on a chain of length 12. The L — 3 fi- 
delity decrease for /3 < — | can be attributed to a one-off effect 
in our renormalization map, as described in Fig. 3. In (b) the 
buffering probability is normalized by the AKLT value 1 /'i L ; 
the L = 9 factor decays to roughly 10 -5 as f3 — ► — 1, though not 
plotted explicitly due to space constraints. Rotation fidelity 
is computed by entangling the chain with a ficticious termi- 
nation spin on the right edge and calculating the overlap of 
the measured chain with the ideal output state. The unique, 
entangled ground state of the doubly-terminated chain can be 
quickly found using sparse matrix methods. This overlap can, 
in turn, be evaluated using the expectation of string opera- 
tors, which for rotation by 6 about the I axis (measurement 
in the {x',y',z} basis), are E x / ® o~ x and E z ® o z . The ini- 
tial state is an eigenstate of S z ® a z and remains so after the 
measurement, implying the square of the fidelity of the output 
state with the ideal state is given by F 2 — |(1 + (E x i <g> a x )). 

Short-ranged variations and buffered logical gates. — 
The reason for the reduced fidelity can be understood 
by qualitatively examining the variations in the ground 
state away from the AKLT point. The left-terminated, 
lcngth-iV AKLT chain has the exact MPS description 

|£(-§)) x E{6 3 } IM ' ' ' l & i> ® O&jv ® • • • ® CfcJ !</>)> for 
\(jf) an arbitrary spin-^ state, o~k the Pauli spin operators 
(except a v =cr x a z ), and bj € {x, y, z} for j £ {1, . . . , N}. 

Observe that the physical swap operation of the spins 
j and j+1 is <Sj*,j+i := Sj • Sj+i + (Sj ■ Sj+1) 2 — 1, so that 
changing the relative weight (3 in Eq. (1) essentially cor- 
responds to a perturbation by swaps. Heuristically, the 
ground state \G(/3)) for j3 rj — | is a coherent superpo- 
sition of the AKLT state and various partially swapped 
AKLT states with reordered matrix products. Standard 
basis measurements are unaffected since different order- 
ings of the anticommuting Pauli operators differ at most 
by a factor —1; this alters the probability of the different 
measurement sequences but not the resulting state. How- 
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ever, consider a non-trivial rotation gate R z (0) resulting 
from measurement outcome \9) on spin j followed by \x) 
on spin The dominant AKLT term yields the log- 

ical action a x (a x R z (9)) — R z {9). But the ordering is 
reversed in the swapped term, and the action on the log- 
ical state is instead (a x R z (9))a x = R z {—6). Interference 
of the two terms reduces the fidelity of the rotation. 

This picture suggests a variation of the measurement 
sequences for logical gates which would mitigate the ef- 
fects of these short-ranged variations (swaps) . Consider 
a block of three spins, in which the R z (9) measurement is 
only attempted on the middle site conditioned on the out- 
comes \z) on the two neighboring sites. The neighboring 
byproducts a z commute with R z {9) and thus, under the 
action of any permutation within the block, the logical 
action is invariant (up to a minus sign). 

This scheme, called buffering, can easily be extended 
to any odd blocklength L. Moreover, buffered measure- 
ments can be concatenated to increase the blocklength; 
Fig 1(c) depicts three L = 3 buffered measurements im- 
plementing an L = 9 measurement. Increasing L yields 
higher gate fidelity, as shown in Fig. 2, though this im- 
provement comes at the expense of the (heralded) success 
probability decreasing exponentially in L, or doubly ex- 
ponentially in the number of concatenated L — 3 steps. 
The buffering overhead is, however, constant in terms of 
the input size of the quantum computation. 

Computational renormalization. — Buffered gate se- 
quences are insensitive to the short-ranged variations ex- 
perienced as one moves away from the AKLT point, but 
still utilize the long-ranged degrees of freedom charac- 
teristic of the Haldane phase. A renormalization group 
can be constructed from these degrees of freedom un- 
der certain coarse-grainings (e.g. [12, 14]), and therefore 
one expects that rotation measurements performed on 
these degrees of freedom, viewed as renormalized spins, 
will have higher rotation fidelity than measurements per- 
formed directly on the physical spins. The challenge in 
MQC is to perform the appropriate renormalized mea- 
surement using only single-site measurements and posts- 
election; for the ID Haldane chain we are guided by the 
heuristic swap analysis. 

Let us elucidate the relation of L=3 buffering to the 
RG. First, we define the map from three spins to one 
which generates the RG flow. Then we explain how 
buffering mimics the desired measurement on the renor- 
malized spin. 

The state of three spin-1 particles can be expressed in 
terms of the total angular momentum of all three spins, 
with J=0, 1,2, 3 components. Due to the symmetry of 
the buffering procedure, it will only take notice of com- 
ponents which are invariant under interchange of the first 
and last spins. This leaves two J=l representations, one 
J=2, and the J=3. It is convenient to think of the re- 
maining J=l sector as a tensor product Hj Cg>%L, where 
Hj carries a J=l irreducible representation of SU(2) and 



is a 2-dimensional 'label' subsystem. An explicit or- 
thonormal basis for this tensor product is 

l°>L - 71 (Ifc)l l*0) 23 + \k) 2 |* )l3 + Ifc> 3 l*o) ia ) , 

|1) L = I (l*>! l*o) 23 " 2 |fc) 2 |* ) 13 + l*>3 l*o) 12 ) , 

where k 6 {x,y,z} and |\&o) = (\ xx ) ~ \w) + \ zz )) is 
the J=0 state of two spin-1 particles. Our RG map is 
then defined as follows: on each sequential block of three 
spin-1 particles, we project the state onto the J=l sector 
and subsequently trace out the label subsystem. 

This yields a new spin-1 state on Wj and, after nor- 
malization, a state on a new spin-1 chain of one-third 
the length. Fig. 3 shows that the map generates an RG 
flow, depicting the state of the block in the Hj <8> Hl 
subspace both before and after one iteration of the map. 
As the system is rotationally-invariant, the state on Hj 
is just the completely mixed state, and we can focus on 
the Bloch vector of the state on H-^. Note that our RG 
scheme does not map the one-parameter family of ground 
states exactly onto itself, but nonetheless flows towards 
states having higher fidelity with the AKLT state. 
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FIG. 3: Bloch vector of the label space "Hi, reduced density 
operator before (dashed) and after (solid) a single RG step 
for —1 < (3 < 1. The state |xs) of Eq. (2) defines the vertical 
axis, while the horizontal corresponds to the superposition 
1+) = Tl(IXs) + for |35>=^(a/5|0) + |1»; the y axis 

is not shown, as the reduced density operator coefficients are 
real in this basis. The Heisenberg antiferromagnet (/3 = 0) is 
denoted HAF. The norm of the Bloch vector, i.e. the radial 
distance from the origin, provides the weight of the symmet- 
ric three-spin J=l sector in the (pre- and post-renormalized) 
ground state \G(/3)), and the projection onto the vertical axis 
indicates the buffering success probability, given successful 
J—l projection. Observe that for /3 > — |, the RG approx- 
imately maps Bloch vectors closer to the AKLT point along 
the original curve parametrized by f3, meaning the correlation 
of the reduced state is effectively renormalized to another /3 
closer to the AKLT point. Meanwhile, the first iteration takes 
/?< — gto/3> — g. Accordingly, we expect that iteration of 
our RG map generates a flow toward the AKLT point. 

It is easy to see that the AKLT state is a fixed point 
of this RG map by using the MPS form. Projecting 
two neighboring spins onto |^o) simply leaves a shorter 
AKLT chain, since a x — —a^ = a"l = 1. Similarly, pro- 
jecting next-neighboring spins onto j^o) again leaves a 
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shorter AKLT chain, with an overall phase of —1. Now 
consider projecting a block of three spins from \Q(— g)) 
onto either |fc)j |0) L or \k)j |1) L . By the above state- 
ments, this yields a shorter AKLT chain, with two spins 
from the block deleted and one projected onto \k). Trac- 
ing out the label subsystem corresponds to mixing these 
two outcomes incoherently. However, this has no effect, 
as the two outcomes are identical. Thus, projecting onto 
the full Hj (g> %l subspace and tracing out the label sub- 
system leaves the form of the AKLT state invariant. 

That 3-spin buffered rotation measurement acts as the 
desired rotation measurement on the renormalized spin 
arises from the following decomposition of the "bare" 
three-spin states in terms of the tensor product Hj ®'Hi J , 

\z, 9, z) 123 oc J\ \6) j |xs) L + J^l component , (2) 
\z,z,z) 123 oc v/| |z)j |0) L + J=^l component, (3) 

where |Xs)— ^g(|0)— v5|l)) is independent of the mea- 
surement angle 9. Buffering implements a projective 
measurement of the label space, rather than a partial 
trace, but because T-Lj and Hl are left unentangled, this 
distinction only affects the success probability and not 
the resulting map. As the weights of the "successful" 
and "failure" outcomes are different, the failure probabil- 
ity of the gate given successful buffering will be different 
from 1/3 away from the AKLT point; interestingly, it can 
actually improve, as shown in Fig. 2(b). In general the 
J=/=l component of the measurement will have non-zero 
overlap with the J^l component of the state \Q(/3)), so 
that the gate fidelity is still less than unity. This is par- 
ticularly relevant for 1, where \G{0)) has increased 
weight on the symmetric three-spin J=2 subspace. 

Conclusion. — We have shown that our renormaliza- 
tion protocol removes the short-ranged variations in the 
(rotationally-invariant) Haldane phase, generating a flow 
toward the AKLT point. Correspondingly, in a practi- 
cal setting our buffering procedure can be used to ensure 
a target gate fidelity (chosen by fault-tolerance consid- 
erations) for any |/3|<1, with attendant decrease in the 
success probability, as shown in Fig. 2. In this sense, the 
quantum computational ability of the spin-chain is a ro- 
bust property of the phase. While we have only given 
results for a single ID chain, it is straightforward to in- 
clude a coupling cphase gate as described in [ ]. Diago- 
nal in the z basis, this gate can thus also be protected by 
buffering. The resulting fidelity improvement almost ex- 
actly follows that of single-qubit rotations as in Fig. 2(a). 

Our quantum computational RG has several unique 
features compared with classical RG methods [ - ]. 
First, it is a renormalization of a class of states rather 
than Hamiltonians [21], though, as Fig. 3 shows, it 
is insensitive to how the "label space" state is treated. 
Remarkably, the map is time-ordered (adaptive in the 
choices of later measurements) , in contrast with the con- 



ventional real-space RG that renormalizes the state or 
Hamiltonian uniformly in space. This is crucial to physi- 
cally implement RG in the present context, both because 
the gates in a quantum circuit provide an implicit time 
ordering, and due to the need to compensate for the in- 
herent randomness of measurement outcomes. 

Focussing on the rotationally invariant Hamiltonians 
of Eq. (1) is well-motivated by physical realizations [20]. 
For example, for spin-1 bosonic atoms trapped in a ID 
optical lattice with tunneling-induced interations, the 
dominant interaction channel is rotationally-invariant s- 
wave scattering. In realizations using microwave induced 
dipole-dipole interacting trapped polar molecules, spher- 
ical symmetry is provided by the choice of polarization, 
intensity, and frequency of the applied fields. 
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